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^ ■ ABSTRACT 

"^ ' Apollonian circle packings arise by repeatedly filling the interstices between mutually tan- 

r^ ■ gent circles with further tangent circles. It is possible for every circle in such a packing to have 

' integer radius of curvature, and we call such a packing an integral Apollonian circle packing. 

I This paper studies number-theoretic properties of the set of integer curvatures appearing in 

^ ■ such packings. Each Descartes quadruple of four tangent circles in the packing gives an inte- 

rs . ger solution to the Descartes equation, which relates the radii of curvature of four mutually 

^^ I tangent circles: x'^ + y'^ + z'^ + w"^ = ^{x + y + z + w)'^. Each integral Apollonian circle packing 

Q>^ ■ is classified by a certain root quadruple of integers that satisfies the Descartes equation, and 

O . that corresponds to a particular quadruple of circles appearing in the packing. We express 

the number of root quadruples with fixed minimal element —n as a class number, and give an 
exact formula for it. We study which integers occur in a given integer packing, and determine 
congruence restrictions which sometimes apply. We present evidence suggesting that the set 
of integer radii of curvatures that appear in an integral Apollonian circle packing has posi- 
tive density, and in fact represents all sufficiently large integers not excluded by congruence 
K^ ■ conditions. Finally, we discuss asymptotic properties of the set of curvatures obtained as the 

^ , packing is recursively constructed from a root quadruple. 
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Apollonian Circle Packings: Number Theory 



1. Introduction 

Place two tangent circles of radius 1/2 inside and tangent to a circle of radius 1. In the two 
resulting curvilinear triangles fit tangent circles as large as possible. Repeat this process for 
the six new curvilinear triangles, and so on. The result is pictured in Figure ^ where each 
circle has been labeled with its curvature — the reciprocal of its radius. 




Figure 1: The integral Apollonian circle packing (-1, 2, 2, 3) 



Remarkably, every circle in Figure ^ has integer curvature. Even more remarkable is that 
if the picture is centered at the origin of the Euclidean plane with the centers of the "2" 
circles on the x-axis, then each circle in the picture has the property that the coordinates of its 
center, multiplied by its curvature, are also integers. In this paper we are concerned with circle 
packings having the first of these properties; the latter property is addressed in a companion 



paper ^| Section 3]. 

An Apollonian circle packing is any packing of circles constructed recursively from an 
initial configuration of four mutually tangent circles by the procedure above. More precisely, 
one starts from a Descartes configuration, which is a set of four mutually tangent circles with 
disjoint interiors, suitably defined. In the example above, the enclosing circle has "interior" 
equal to its exterior, and its curvature is given a negative sign. Recall that in a quadruple of 
mutually touching circles the curvatures (a, b, c, d) satisfy the Descartes equation 

a^ + b^ + c^ + (f = -{a + b + c + df, (1.1) 

as observed by Descartes in 1643 (in an equivalent form). Any quadruple (a, b, c, d) satisfy- 
ing this equation is called a Descartes quadruple. An integral Apollonian circle packing is an 
Apollonian circle packing in which every circle has an integer curvature. The starting point of 
this paper is the observation that if an initial Descartes configuration has all integral curva- 
tures, then the whole packing is integral, and conversely. This integrality property of packings 
has been discovered repeatedly; perhaps the first observation of it is in the 1937 note of F. 
Soddy ^3 "The bowl of integers and the Hexlet" . It is discussed in some detail in Aharonov 
and Stephenson p. 

In this paper we study integral Apollonian circle packings viewed as equivalent under 
Euclidean motions, an operation which preserves the curvatures of all circles. Such packings 
are classified by their root quadruple, a notion defined in §3. This is the "smallest" quadruple 
in the packing as measured in terms of curvatures of the circles. In the packing above the 
root quadruple is (—1, 2, 2, 3), where —1 represents the (negative) curvature of the bounding 
circle. We study the set of integers (curvatures) represented by a packing using the Apollonian 
group A, which is a subgroup of GL{4:,'L) which acts on integer Descartes quadruples. This 
action permits one to "walk around" on a fixed Apollonian packing, moving from one Descartes 
quadruple to any other quadruple in the same packing. The Apollonian group was introduced 
by Hirst [22 in 1967, who used it bounding the Hausdorff dimension of the residual set of 
an Apollonian packing; it was also used in Soderberg,,46^ and Aharonov and Stephenson ^. 
Descartes quadruples associated to different root quadruples cannot be reached by the action 
of A, and the action of the Apollonian group partitions the set of integer Descartes quadruples 
into infinitely many equivalence classes (according to which integral Apollonian packing they 



belong.) By scaling an integer Apollonian packing by an appropriate homothety, one may 
obtain a primitive integral Apollonian packing, which is one whose Descartes quadruples have 
integer curvatures with greatest common divisor 1. Thus the study of integral Apollonian 
packings essentially reduces to the study of primitive packings. 

The simplest integral Apollonian circle packing is the one with root quadruple (0, 0, 1, 1), 
which is pictured in Figure 2. This packing is special in several ways. It is degenerate in 
that it has two circles with "center at infinity", whose boundaries are straight lines, and it is 
the only primitive integral Apollonian circle packing that is unbounded. It is also the only 
primitive integral Apollonian circle packing that contains infinitely many copies of the root 
quadruple. This particular packing has already played a role in number theory. That part of 
the packing in an interval of length two between the tangencies of two adjacent circles of radius 
one, consisting of the (infinite) set of circles tangent to one of the straight lines, forms a set of 
"Ford circles", after shrinking all circles by a factor of two. These circles, introduced by Ford 
(see ^ni)) can be labelled by the Farey fractions on the interval (0, 1) and used to prove basic 
results in one-dimensional Diophantine approximation connected with the Markoff spectrum, 
see Rademacher |SH1 PP- 41-46] and Nicholls PSJ . 

In this paper our interest is in those properties of the set of integral Apollonian circle 
packings that are of a Diophantine nature. These include the distribution of integer Descartes 
quadruples, of integer root quadruples, and the representation and the distribution of the 
integers (curvatures) occurring in a fixed integral Apollonian circle packing. Finally we consider 
the size distribution of elements in the Apollonian group, a group of integer matrices associated 
to such packings. 

To begin with, the full set of all integer Descartes quadruples (taken over all integral 
Apollonian packings) is enumerated by the integer solutions to the Descartes equation (up 
to a sign.) In §2 we determine asymptotics for the total number of integer solutions to the 
Descartes equation of Euclidean norm below a given bound. 

In §3 we define the Apollonian group. We describe a reduction theory which multiplies 
Descartes quadruples by elements of this group and uses it to find a quadruple of smallest size 
in a given packing, called a root quadruple. We prove the existence and uniqueness of a root 
quadruple associated to each integral Apollonian packing. 




Figure 2: The Apollonian circle packing (0, 0, 1, 1). 



In §4 we study the root quadruples of primitive integer packings. We show that the number 
N*g^^{—n) of primitive root quadruples containing a given negative integer —n is equal to 
the class number of primitive integral binary quadratic forms of discriminant — 4n^, under 
GL(2,Z)-equivalence. Using this fact we obtain in Theorem 14.31 an exact formula for the 
number of such quadruples, which is 

Kooti-n) = 5 11(1 - ^^) + 2-(")-^-\ (1.2) 

p\n 

where X-4('^) = (— 1)("'~^)/2 for n and for even n, a;(n) denotes the number of distinct primes 
dividing n, and (5„ = 1 if n = 2( mod 4) and 5„ = otherwise. This formula was independently 
discovered and proved by S. Northshield [HZ]. As a consequence we show that N*^^^{—n) has 
an upper bound 0(n log log n) and a lower bound ^( nogiogn) )' respectively. These bounds are 
sharp up to multiplicative constants. 

In §5 we study the integer curvatures appearing in a single integral Apollonian packing, 
counting integers with multiplicity. D. Boyd [Zj showed that the number of circles occurring in 
a bounded Apollonian packing having curvature less than a bound T grows like T°''^°^^' , where 
a ~ 1.30... is the Hausdorff dimension of the residual set of any Apollonian circle packing. 
This result applies to integral Apollonian packings. We observe that these integers can be put 
in one-to-one correspondence with elements of the Apollonian group, using the root quadruple. 
Using this result we show that the number of elements of the Apollonian group which have 
norm less than T is of order T'^+^W ^ ^s T — > oo. 

In §6 we study the integer curvatures appearing in a packing, counted without multiplicity. 
We show that there are always nontrivial congruence restrictions {mod 24) on the integers that 
occur. We give some evidence suggesting that such congruence restrictions can only involve 
powers of the primes 2 and 3. We conjecture that in any integral Apollonian packing, all 
sufficiently large integers occur, provided they are not excluded by a congruence condition. This 
may be a hard problem, however, since we show that it is analogous to Zaremba's conjecture 
stating that there is a fixed integer K such that for all denominators n>2 there is a rational 
- in lowest terms whose continued fraction expansion has all partial quotients bounded by K. 

In §7 we study the set of integer curvatures at "depth n" in an integral Apollonian packing, 
where n measures the distance to the root quadruple. There are exactly 4x 3"^^ such elements. 



We determine the maximal and minimal curvature in this set, and also formulate a conjecture 
concerning the asymptotic size of the median curvature as n ^ cxd. These problems are related 
to the joint spectral radius of the matrix generators S = {S*!, 5*2, §3, 5*4} of the Apollonian 
group, which we determine. 

In §8 we conclude with some directions for further work and a list of open problems. 

There has been extensive previous work on various aspects of Apollonian circle packings, 
related to geometry, group theory and fractals. The name "Apollonian packing" traces back 
to Kasner and Supnick [231 ^^ 1943. It has been popularized by Mandlebrot [SHI P- 169ff], who 
observed a connection with work of Apollonius of Perga, around 200BC, see also Coxeter jH]. 
Further discussion and references can be found in Aharonov and Stephenson 1_ and Wilker |S2] • 
See also the companion papers ^7\, JH1)CH1 and [H], where we investigate a variety of group- 
theoretic properties of these configurations, as well as various extensions to higher dimensions 
and other spaces, such as hyperbolic space. 

Notation. We use several different measures of the size of a vector w = {wi, 102,103,104). 
In particular we set L(w) := wi + 102 + W3 + w^, a measure of size which is not a norm. 
We let H{-w) := {wi + W2 + w"^ + w^y'^ or |w| denote the Euclidean norm, while ||w||oo := 
maxi<j<4 \wi\ denotes the supremum norm of w. 

2. Integral Descartes Quadruples 

An Apollonian circle packing is integral if every circle of the packing has an integer curvature. 
From (|1.1() it follows that if a,b, c, are given, the curvatures d, d! of the two circles that are 
tangent to all three satisfy 

d, d' = a + b + c^ '^qabci 

where 

qabc = Vab + bc + ac. 

Hence 

d + d' = 2{a + b + c). (2.1) 



In other words, given four mutually tangent circles with curvatures a, b, c, d, the curvature of 
the other circle that touches the first three is given by 

d' = 2a + 2b + 2c- d. (2.2) 

It follows that an Apollonian packing is integral if the starting Descartes quadruple consists 
entirely of integers. 

The relation H2.1() is the basis of the integrality property of Apollonian packings. It gener- 
alizes to n dimensions, where the curvatures Xi of a set of n + 1 mutually tangent spheres in 
M" (having distinct tangents) are related to the curvatures Xq and Xn+2 of the two spheres 
that are tangent of all of these by 

2 

-'^O + Xn+2 = ri^l + X2 + ■ ■ ■ + Xn). 

n — 1 
This relation gives integrality in dimensions n = 2 and n = 3; the three dimensional case 
is studied in Boyd :5 . It even generalizes further to sets of equally inclined spheres with 
inclination parameter 7, with the constant — ^; the case 7 = — 1 is the mutually tangent case, 
cf. Mauldon |3T] and Weiss ^ Theorem 3]. 

Definition 2.1. (i) An integer Descartes quadruple a = {a,b,c,d) G Z^ is any integer repre- 
sentation of zero by the indefinite integral quaternary quadratic form, 

Qviw, X, y, z) := 2{w'^ + x"^ + y"^ + z"^) - {w + x + y + z)"^, 

which we call the Descartes integral form. That is, writing v = {w, x, y, z) , we have 
Qv{w, X, y, z) = v^Qx'V, where 



Q 



V 



1 -1 -1 -1 

-1 1 -1 -1 

-1 -1 1 -1 

-1 -1 -1 1 



(2.3) 



This quadratic form has determinant det{Qx>) = — 16 and, on identifying the form Qj) with its 
symmetric integral matrix, it satisfies Q^, = 4/. 

(ii) An integer Descartes quadruple is primitive if gcd{a,b,c,d) = 1. 

In studying the geometry of Apollonian packings (JZI-jn]) |21]) we use instead a scaled version 
of the Descartes integral form, namely the Descartes quadratic form Q2 '■= 2Qt>- 



Definition 2.2. The size of any real quadruple (a, b, c, d) G M^ is measured by the Euclidean 
height H(a), which is: 



i/(a):=(a2 + 62 + 0^ + ^2)1/2. 



(2.4) 



Now let Nx){T) count the number of integer Descartes quadruples with Euclidean height 
at most T. We shall relate this quantity to the number Nc{T) of integer Lorentz quadruples 
of height at most T, where Lorentz quadruples are those quadruples that satisfy the Lorentz 
equation 



W^ + X^ + Y^ + Z^ = 0. 



These are the zero vectors of the Lorentz quadratic form 



(2.5) 



Qc{W, X, Y, Z) = -W'^ + ^2 + y2 + ^2^ 



(2.6) 



whose matrix representation is 



Qn 



-10 
10 
10 
1 



Similarly we shall relate the number of primitive integer Descartes quadruples, denoted N^{T), 
to the number of primitive integer Lorentz quadruples of height at most T, denoted N^{T). 
We show that there is a one-to-one height preserving correspondence between integer Descartes 
quadruples and integer Lorentzian quadruples. Introduce the matrix Jq defined by 
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1 


1 


1 


1 


-1 


-1 


-1 


1 


-1 


-1 


-1 


1 



(2.7) 



and note that Jq = /. The Descartes and Lorentz forms are related by 



Qt> = 2Jq QcJq, 



(2.8) 



which leads to a relation between their zero vectors. 



Lemma 2.1. The mapping (W, X, Y, Z) = Jo{w, x, y, z) gives a bijection from the set 
{w, X, y, z) of real Descartes quadruples to that of real Lorentz quadruples {W, X, Y, Z) 
which preserves height. It restricts to a bijection from the set of integer Descartes quadruples 
to integer Lorentz quadruples, so that Nx>{T) = Nc{T), for all T > 0, and from, primitive 
integer Descartes quadruples to primitive integer Lorentz quadruples, so that N^{T) = N^(T), 
for all T > 0. 

Proof. An easy calculation shows that the mapping takes real solutions of one equation to 
solutions of the other and that the inverse mapping is (tt;, x, y, z)'^ = Jo{W, X, Y, Z)^, so that 
it is a bijection. The mapping takes integer Descartes quadruples to integer Lorentz quadruples 
because any integer solution to the Descartes equation satisfies w + x + y + z = (mod 2). This 
also holds in the reverse direction because integer solutions to the Lorentz form also satisfy 
W + X + Y + Z = {mod 2), as follows by reducing ()2.5|) {mod 2). It is easy to check that 
primitive integer Descartes quadruples correspond to primitive integer Lorentz quadruples. □ 
Counting the number of integer Descartes quadruples of height below a given bound T is 
the same as counting integer Lorentz quadruples. This is a special case of the classical problem 
of estimating the number of representations of a fixed integer by a fixed diagonal quadratic 
form, on which there is an enormous literature. For example Ratcliffe and Tschantz |4nj 
give asymptotic estimates with good error terms for the number of solutions for the equation 
j^2 _j_ y2 _|_ _^2 _ ^2 _ ^^ q£ Euclidean height below a given bound, for all k ^ 0. (They treat 
Lorentzian forms in n variables.) Rather surprisingly the case k = seems not to be treated 
in the published literature, so we derive an asymptotic formula with error term below. The 
main term in this asymptotic formula was found in 1993 by W. Duke jHI (unpublished) in the 
course of establishing an equidistribution result for its solutions. 

Theorem 2.2. The number of integer Descartes quadruples Nx>{T) of Euclidean height at 
most T satisfies Nt>(T) = Nc(T), and 

^^(^) = JTTV T ^' + 0{T{logTf), (2.9) 



as T —> oo, in which 

L(2, x-4) = y T^ — -T^ ~ 0. 



9159. 

Vln + 1 V 
n.=0 



The number N^ (T) of primitive integer Apollonian quadruples of Euclidean height less than T 
satisfies N^{T) = Nl{T) and 

Nl{T) = ^ T' + 0{T{lo9Tf), (2.10) 

as T —> oo. 

Proof. By Lemma l2. II it suffices to estimate Nc{T). Let ^3(771) denote the number of integer 
representations of ttt, as a sum of three squares, ahowing positive, negative and zero integers. 
Rewriting the Lorentz equation as X'^ + Y^ + Z^ = W'^ we obtain for integer T that 

T 

Nc{V2T) = l + 2Y,r3{m^), (2.11) 

m=l 

since there are two choices for W whenever W ^ 0. A general form for r-^^m) was obtained in 
1801 by Gauss |161 Articles 291-292], while in the special case r^^m?) a simpler form holds, 
given in 1906 by Hurwitz 22 . This is reformulated in Sandham j43l p. 231] in the form: if 
m = YipP'^''^"^^ ^-iid p runs over the primes and rriodd = m2~'^'^^"^> , then 

-I— r 2 2 

rsim^) = Qruodd JJ (H h .... + 



p rf-pV"-) 

p=3 (mod 4) ^ 

« n {- TT^ -]■ P-12) 

p=l (mod 2) \ / 



Sandham observes that this formula is equivalent to 



F{s) := y '-^^ = 6(1 - 2i-)^M(^ (2.13) 



where 



00 . QO 



L{s, x-4) := E (— )m- = Y, ^ ' (2.14) 

The right hand side of (|2.13j) is a meromorphic function in the s-plane, which has a simple 
pole at s = 2 with residue 

3C(2) TT^ 

ci - - 



L(2,x-4) 2L(2,x-4)' 
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and has no other poles for Ks > 1. One could now proceed by a standard contour integral 
approach ^ to obtain Nc{\2,T) equals ciT^ plus a lower-order error term, so that Nc{T) = 
gCiT^ + o(T'^). However we will derive (|2.9|) directly from the exact formula p.l2|) . We have, 

Nc{V2T) = 1 + 2 Yl E ^3((2n-l)2) 

'-2J^2J ^ 



1+12 E E (2r^-i) n (1+-+--+ 



i<j<log2T n=l P=3 (mod 4) 

p|2n-l 



p pep(2n-l)- 



Expanding the products above and using X]n=i(2^ — 1) = f/^ + 0{U), one obtains 



NciV2T) = 1 + 12EE7^IE E (2"^-im 

k>0 Pk '^ \j>l "1=1 

/ \ 

12 



4 



E^^ E Wp^r' + of(E2'=EK)^V ('•'') 

fc>0 J>o.^fc ^ / \ fc>0 Pfc<T *= / 



\ Pfc<T/2J / 



in which Pk = p^ ...p^j^ with all pi = 2> [mod 4) and all ei > 1. If the condition P^ < T/2^ 
were dropped in the first sum in parentheses above, then it would sum to X72^^' ^^ '-'^^ ^^^^ 
by examining the associated Euler product, which is 

^ "^ ^'^ p=3 (mod 4) ^ m=l 

evaluated at s = 2, since jzr^rj = 1 + 2p^^ + 2p^^* + ... The error introduced by truncating 
this Dirichlet series at s = 2 at n < T is bounded by 



-a(m)^|^2-M^-d(^^ ^ logT 



m=T m=T m=T 

in which u{n) counts the number of distinct prime divisors of m (without multiplicity), d{n) 
counts the number of divisors on n, and the last estimate uses ^ the fact f j2Ul Theorem 
319]) that the average order of d{n) is logn. The remainder term in (|2.15|) is bounded by 

^Integrate F{s) against —ds on a vertical contour 'Rs = c for some c with 1 < c < 2, cf. Davenport |12l 
Chapter 17]. 

*Sum in blocks 2-'T < n < 2-'^^T, viewing denominators as nearly constant and averaging over numerators. 
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0(T(logT)2), because 



A:>0 Pk<T m=l m=l m=l 



Since C(2) = ^, these estimates combine to give 

Nc{V2T) = f T' + 0{T {logTf). 

2L{2,x-4) 

Rescaling T yields ()2.9() . 

To handle primitive Lorentz quadruples, we use the function r^{m) which counts the number 
of primitive integer representations of ?7i as a sum of three squares, using positive, negative 
and zero integers. Then r3(m^) = ^^ i ^r^id?), which yields 

The function F*{s) is analytic in 3f?(s) > 1 except for a simple pole at s = 2 with residue -^-ci, 
which gives the constant in the main term of (|2.1()|) . To obtain the error estimate (|2.1()j) one 
can use (|2.9|) and Mobius inversion, with r3(?TT,^) = J2d I rra/^('^)^3((T)^)- ^^ omit details. □ 

Remarks. (1) Various Dirichlet series associated to zero solutions of indefinite quadratic 
forms have meromorphic continuations to C, cf. Andrianov 2j. These can used to obtain 
asymptotics for the number of solutions satisfying various side conditions. 

(2) The real solutions of the homogeneous equation Qc{w, x, y, z) = —w'^ + x^ + y^ + z^ = 
form the light cone in special relativity. 

3. Reduction Theory and Root Quadruples 

In this section we describe, given an Apollonian circle packing and a Descartes quadruple in 
it, a reduction procedure which, if it halts, identifies within it a unique Descartes quadruple 
(a, h, c, d) which is "minimal" . This quadruple is called the root quadruple of the packing. This 
procedure always halts for integral Apollonian packings. 

Definition 3.1. The Apollonian group A is the group generated by the four integer 4x4 



12 



matrices 



Si 



-12 2 2 
10 
10 
1 



So 



10 
2-122 
10 
1 



S. 



10 
10 
2 2-12 
1 



04 



10 
10 
10 
2 2 2-1 



As mentioned earlier, the Apollonian group was introduced in the 1967 paper of Hirst PT] . 
and was later used in Soderberg |1S] and Aharonov and Stephenson ^ in studying Apollonian 
packings. 

We view real Descartes quadruples v = (a, b, c, d)'^ as column vectors, and the Apol- 
lonian group A acts by matrix multiplication, sending v to Mv, for M £ A. The action 
takes Descartes quadruples to Descartes quadruples, because A C Autz{Qv), the set of real 
automorphs of the Descartes integral quadratic form Qx> given in H2.3|) . That is, each such M 
satisfies 

M^QvM = Qv, for ah M e A, 

a relation which it suffices to check on the four generators Si £ A. 

The elements Sj have a geometric meaning as corresponding to inversion in one of the 
four circles of a Descartes quadruple to give a new quadruple in the same circle packing, as 
explained in \n\ Section 2]. That paper showed that this group with the given generators is a 
Coxeter group whose only relations are Sf = S2 = S^ = S^ = I. 

The reduction procedure attempts to reduce the size of the elements in a Descartes quadru- 
ple by applying one of the generators Sj to take the quadruple v = (a, b, c, d) viewed as a 
column vector to the new quadruple SjV, until further decrease is not possible. Each packing 
has a well-defined sign, which is the sign of L{v) := a + b + c + dof any Descartes quadru- 
ple in the packing. (The sign is independent of the quadruple chosen in the packing; see 
Lemma l3.ir ii) below.) We describe the reduction procedure for packings with positive sign, 
those with L{v) := a + b + c + d > 0; the reduction procedure for negative sign packings is 
obtained by conjugating by the inversion {a,b,c,d) — > {—a,—b,—c,—d), which takes negative 
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sign packings to positive sign packings. Thus we suppose L{v) ■.= a + b + c + d>0 and for 
simplicity consider the case where the quadruple is ordered a <b < c < d. We consider which 
SjV can decrease the sum L{v) = a + b -\- c + d, and show below this can never occur using 
Si, S2 or 53. Note that ^4(0, b, c, d)^ = (a, b, c, d')^ where d' = 2{a + b + c) — d. 

Lemma 3.1. Suppose that v = (a, b, c, d) is a real Descartes quadruple. Let v have its 
elements ordered a < b < c < d, and set d' = 2{a + b + c) — d. 

(i) If L(v) = a + b+ c+ d > 0, then a + b >0, with equality holding only if a = b = and 
c = d. As a consequence, we always have b > 0. 

(ii) If L{v) = a+b+c+d > 0, then each L{Sjv) > 0. In particular L{S 4V) = a+b+c+d' > 0. 

(Hi) If a > 0, so that L{v) = a + b+c+d>0, then d' < c < d. If d' < c then the matrix 
S4 that changes d to d' strictly decreases the sum L{v) := a + b + c + d, and it is the only 
generator of A that does so. If d' = c then necessarily c = d = d' and no generator Si of A 
decreases L{v). 

Proof, (i) If a > then we are done, so assume a < 0. Suppose first that < b < c < d. 
Let X = —{a + b), y = —ab, so that y > 0. In terms of these variables, the Descartes equation 
becomes 

2(x2 + 2y + c^ + d^) = {-x + c + df, 

which simplifies to 

[c-df + x^ + 'iy + 2{c + d)x = Q. (3.1) 

This equation cannot hold if 2; > 0, because all terms on the left are nonnegative, and some 
are positive, hence a + 6 > 0. If x = 0, then (|3.1() implies y = and c = d, whence a = b = 
and c = d. 

The remaining case is a < 6 < 0. We first observe that in any Descartes quadruple, at least 
three terms have the same sign. Indeed the Descartes equation can be rearranged as 

{a-bf + {c-d)^ = ia + b)ic + d), 

and iia<b<0<c<d then the left side is nonnegative while the right side is strictly 
negative, a contradiction. Thus we must have a<b<c<0<d, since a + b + c + d > 0. 
Now the Descartes quadruple {—d, — c, —b, —a) has —d < < — c, so the preceding case applies 
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to show that —d — c > 0. We conclude that a + b + c + d < 0, contradicting the fact that 
a + b + c + d > 0. This completes the proof of (i). 

(ii) The ordering of entries gives L(S'iv) > L(52v) > L(53v) > L(S'4v), so it suffices to 
prove L(54v) > 0. The Descartes equation implies that 



d,d' = a + b + ci: 2qabc, where qabc = Va6 + bc + 



ca. 



We have a + b + c + d' = 2(a + 6 + c) — 2\/a& + bc + ac > because a + b + c> (using (i)) and 

{a + b + cf - {ab + bc + ac) = -{{a + bf + (6 + cf + (a + cf) > 0. 

(iii) The Descartes equation (|1.1|) gives 

d' = a + b + c— 2\^ab + bc + ac. 

Thus 

d'-c = a + b- 2yjab + bc + ac <a + b- i/4(a + b)c <a + b- ^/{a + b)'^ = 0. 

If d' < c < d then the sum L{v') = a + b -\- c + d' < L{v), so the sum decreases. If Si 
changes c to c' , then c' = 2(a + b + d) — c> 2{a + b + c) — c> c because a + 6 > by (i), so the 
sum L{v) does not decrease in this case. Similarly the sum does not decrease if Si changes b 
to b' or a to a'. In the case of equality d' = c, one easily checks that c = d = d' , which forces 
a = b = 0, and no Si decreases the sum -L(v). □ 

Definition 3.2. A Descartes quadruple v = (a, b, c, d) with L{v) =a + 6 + c + (i>0isa root 
quadruple ifa<0<b<c<d and a + b + c> d. 

Note that the last inequality above is equivalent to the condition d' = 2{a + b + c) — d > d. 

Reduction algorithm. 

Input: A real Descartes quadruple (a, b, c, d) with a + b + c + d > 0. 

(1) Test in order 1 < i < 4 whether some Si decreases the sum a + b + c + d. It so, apply 
it to produce a new quadruple and continue. 
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(2) If no Si decreases the sum, order the elements of the quadruple in increasing order and 
halt. 

The reduction algorithm takes real quadruples as input, and is not always guaranteed to 
halt. The following theorem shows that when the algorithm is given an integer Descartes 
quadruple as input, it always halts, and outputs a root quadruple. In the algorithm, the 
element Si that decreases the sum necessarily decreases the largest element in the quadruple, 
leaving the other three elements unchanged. The proof below establishes that in all cases where 
a reduction is possible, the largest element of the quadruple is unique, so that the choice of Si 
in the reduction step is unique. There do exist quadruples with a tie in the largest element, 
such as (0, 0, 1, 1), but the vector (a, 5, c, d) cannot then be further reduced. 

Theorem 3.2. (1) // the reduction algorithm ever encounters some element a < 0, then it 
will halt at a root quadruple in finitely many more steps. 

(2) // o, 6, c, d are integers, then the reduction algorithm will halt at a root quadruple in 
finitely many steps. 

(3) A root quadruple is unique if it exists. However an Apollonian circle packing may 
contain more than one Descartes configuration yielding this quadruple. 

Proof. 

(1) Geometrically a Descartes quadruple with a < describes a circle of of radius 1/a 

enclosing three mutually tangent circles of radii 1/6, 1/c, \/d. All circles in the packing lie 

inside this bounding circle of radius 1/a. Each non-halting reduction produces a new circle of 

radius 1/d' > 1/d, which covers an area of vr/d' , and this is at least vr/d^. Since there is a 

total area of vr/a^ which can be covered, and all circles except the one with radius 1/a have 

disjoint interiors, this process must halt in at most (— ) steps. 

L a J 

(2) Let qabc = \/a6 + bc + ac = {a + b + c — d)/2 G N. After each reduction, the sum 
a + h + c + d decreases by 4gafec- By Lemma IXTl the sum a + h-\- c + d\s bounded below by 0. 
Therefore this process halts after finitely many steps. 

(3) If (a, h, c, d) is a root quadruple of an Apollonian packing, then the numbers a, 6, c, d 
are the curvatures of the largest circles contained in this packing, hence they are unique. On 
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the other hand, the Apollonian packing may contain more than one copy of this quadruple, for 
example, (—1,2,2,3) appears twice in the packing shown in Figure ^ and (0,0,1,1) appears 
infinitely many times in the packing in Figure ^ □ 

Root quadruples lead to a partition of the set (5(^) of all integer Descartes quadruples. 
This set partitions into Q{Z)+ U {(0,0,0,0)} U Q{I.)-, where 

Q(Z)+ = {(a, b, c, d) £ Q{Z) -.a + b + c + dyO} (3.2) 

and Q(Z)^ = — (5(Z)+. Next we have the partition 

oo 

Q(Z)+ = U kQ{Z)+,^ , (3.3) 

fc=i 

where Q{'Z)prim enumerates all primitive integer Descartes quadruples in Q(Z)"'". These latter 

are exactly the Descartes quadruples occurring in all primitive integer Apollonian packings, so 

we may further partition Q(Z)i.^^ into a union of the sets Q{Vv), where Q{Pv) denotes the 

set of all Descartes quadruples in the circle packing V^i having primitive root quadruple V, i.e. 

Q{^)Um= U ^(^^)- (3.4) 

primitive root 
quadruple V 

We study the distribution of root quadruples in §4 and the set of integers in a given packing 
Vt) in §5 and §6. 

By definition the Apollonian group labels all the (unordered) Descartes quadruples in a 
fixed Apollonian packing. We now show that it has the additional property that for a given 
integral Apollonian packing, the integer curvatures of all circles not in the root quadruple lie 
in one-to-one correspondence with the non-identity elements of the Apollonian group. 

Theorem 3.3. Let Vv be the integer Apollonian circle packing with root quadruple v = 
(a, b, c, d) . Then the set of integer curvatures occurring in V, counted with multiplicity, 
consists of the four elements of v plus the largest elements of each vector Mv, where M runs 
over all elements of the Apollonian group A. 

Proof. Let M = Si^ • • • Sj^ be a reduced word in the generators of A, that is Sj^, ^ 5'jj.^^ for 
1 < k < n. The main point of the proof is that if -w*-"' = S-i^ ■ ■ ■ Si-^\r, then w'"-* is obtained 
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from w'""^) by changing one entry, and the new entry inserted is always the largest entry in 
the new vector. (It may be tied for largest value.) We prove this by induction on n. In the base 
case n = 1, there are four possible vectors SiV, whose inserted entries are a' = 2(6 + c + d) — a, 
b' = 2{a + c + d) — b, c' = 2{a + b + d) — c, and d' = 2{a + b + c) — d, respectively. Since 
a < b < c < d we have d' < c' < b' < a' , and since v is a root quadruple, we have d' > d, as 
asserted. 

For the induction step, where n >2, there are only three choices for Si„ since Si^ 7^ 'S'i„_i • If 
the elements of w*-"""'^^ are labelled in increasing order as Wi < W2 < 1^3 < W4 , 
then we may choose the labels (in case of a tie for the largest element) so that w\ was 
produced at step n — 1, by the induction hypothesis. Thus exchanging u;^ is forbidden at 

(n) 

step n, hence if w\ denotes the new value produced at the next step, then 

4"^ > 2(u;S"-^) + 4"-^) + u;("-^) - u;J"-') 

> 2u;!"-') + 24"-^) + 4"-^) > u;i"-'\ (3.5) 

because w^ + W2 > by Lemma IXlT i). This completes the induction step. 

The inversion operation produces one new circle in the packing, namely the new value 
added in the Descartes quadruple, and ()3.5|) shows that its curvature is ||Mv||oo, where ||w||oo 
denotes the supremum norm of the vector w. 

Every circle in the packing is produced in this procedure, by definition of the Apollonian 
group. That all non-identity words M £ A label distinct circles is clear geometrically from the 
tree structure of the packing. □ 

4. Distribution of Primitive Integer Root Quadruples 

In this section we count integer Apollonian circle packings in terms of the size of their root 
quadruples. Recall that a Descartes quadruple (a, b, c, d) is a root quadruple if a < < 6 < 
c < d and d' = 2{a + b + c) — d > d > 0. It suffices to study primitive packings, i.e. ones whose 
integer quadruples are relatively prime. We begin with a Diophantine characterization of root 
quadruples. 

Theorem 4.1. Given a solution {a,b,c,d) G Z^ to the Descartes equation 

{a + b + c + d)'^ = 2(0^ + b'^ + c^ + d'^), 
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define {x,di,d2,m) by 



a 




"10 


" 




X 




X 




b 




-1 1 







di 




di — X 




c 




-1 1 







d2 




d2 - X 




d _ 




-111 


—2 




m 




-2m + di + d2- 


- X 



(4.1) 



Then {x,di,d2,m) G Z^ satisfies 



X +m = did2 



(4.2) 



Conversely, any solution {x,di,d2,Tn) £ Z'* to this equation yields an integer solution to the 
Descartes equation as above. In addition: 

(i) The solution {a,b,c,d) is primitive if and only if gcd{x,di,d2) = 1. 

(a) The solution (a, b, c, d) with a + b + c + d>Oisa root quadruple if and only if 
X <0 <2m < di < d2. 



Proof. The first part of the theorem requires, to have m € Z, that a + b + c + d = (mod 2). 
This follows from the Descartes equation by reduction (mod 2). 

For (i), note that gcd(x, di, ^2) = gcd(a, 6, c) = gcd{a,b,c,d). 

For (ii), the condition a<0<b<c<d implies successively x < 0, di < d2, di — 2x = 
b—a > 0, and —2m+di = d—c > 0. Finally the root quadruple condition d' = 2{a+b+c) — d > 
d > gives d' = 2m > 0. Thus x < < 2m < di < d2. The converse implication follows 
similarly. □ 

We now study primitive integer root quadruples with a = —n, for n E Z>o. Theorem 14. II 
shows that they are in one-to-one correspondence with the integer solutions {m,di,d2) to 



n + m = did2, 



(4.3) 



< 2m < di < d2 and gcd(n, di, ^2) = 1 . 



(4.4) 



For each of n = 0, 1, 2, there is one primitive root quadruple with a = —n, namely, (0, 0, 1, 1), 
(—1,2,2,3), (—2, 3, 6, 7), respectively. For n = 3, there are two, (—3,4,12,13) and (—3, 5, 8, 8). 
As an example of a nonsymmetric integral Apollonian circle packing. Figure 3 pictures the 
packing (—6, 11, 14, 15). 
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Figure 3: The nonsymmetric packing (—6, 11, 14, 15). 
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Let N*ggf.{—n) denote the number of primitive root quadruples with negative element —n. 
One has N*g^^{—n) > 1 for all n > 0, since (x, di, d2, m) = (— n, l,n^,0) in Theorem 14.11 
produces the primitive root quadruple (a, h, c, d) = (— n, n + 1, n{n + 1), n(n + 1) + 1). Table 
1 and Table 2 present selected values of N*^^f.{—n). 
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Table 1: N*^^^{-n) for small n 



p 


N{-p) 


P 


N{-p) 


P 


N{-p) 


P 


N{-p) 


1009 


253 


3001 


751 


4007 


1003 


5011 


1254 


1013 


254 


3011 


754 


4013 


1004 


10007 


2503 


2003 


502 


4001 


1001 


5003 


1252 


10009 


2503 


2011 


504 


4003 


1002 


5009 


1253 


20011 


5004 



Table 2: N*^^^{—p) for selected prime p. 



Examination of numerical data led the the third author and S. Northshield independently 



to conjecture an exact formula for N*^ 



rooty 



-n] 



stated as Theorem 14.31 below. Northshield [37] 



obtained a proof, and we give another here. We establish this result by showing that N*.^^i{—n) 
can be interpreted as a class number, namely N*^^^{—n) = h^{—An?) where for — 4A < the 
quantity /i^(— 4A) counts the number of GL ( 2, Z) -equivalence classes of primitive integral 
binary quadratic forms 

[A, B, C] := AT^ + 2BTU + CC/^ 

of discriminant — 4A = AB^ — AAC; a form is primitive if gcd{A,B,C) = 1. By comparing 
GL(2,Z)-actions and 5'L(2,Z)-actions on positive definite binary quadratic forms one obtains 
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the relation 

/i±(-4A) = ^(/i(-4A) + a(-4A)) (4.5) 

where /i(— 4A) is the usual binary quadratic form class number (for SL(2, Z)-equivalence) and 
a(— 4A) counts the 5L ( 2, Z) -equivalence classes containing an ambiguous reduced binary form 
of discriminant — 4A (defined below). The exact formula in Theorem 14.31 is deduced from 
classical formulas for /i(— 4n^) and a(— 4n^). 

Theorem 4.2. The Apollonian root quadruples {—n, x, y, z) with —n <0<x<y<z are in 
one-to-one correspondence with positive definite integral binary quadratic forms of discriminant 
—All? having nonnegative middle coefficient. The associated reduced binary quadratic form 
[A, B, C] = AT^ + 2BTU + CU^ is given by 

[A,B,C] :=[-n + x, -(-n + x + y - z), -n + y]. 

Primitive root quadruples correspond to reduced binary quadratic forms having a nonnegative 
middle coefficient. In particular, the number N*^^^{—n) of primitive root quadruples with least 
element —n satisfiies 

where h (—An^) is the number of GL{2,'Z)- equivalence classes of positive definite primitive 
binary integral forms of discriminant —4n^. 

Remark. A positive definite form [^,i3,C] is reduced if < \2B\ < A < C. Reduced forms 
with nonnegative middle coefficients enumerate form clases under GL{2, Z)-action, rather than 
the S'L(2,Z)-action studied by Gauss; the GL(2,Z)-action corresponds to the notion of equiv- 
alence of (definite) quadratic forms used by Legendre. 

Proof. Recall that the Descartes quadratic form is 

Qv{w, x,y,z) := {w -\-x -\-y -\- z)^ - 2(w^ + x^ + y^ + z^). 
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A Descartes quadruple is any integer solution to QT>iw, x, y, z) = and a Descartes quadruple 
is primitive if gcd{w, x,y,z) = 1. Any integer solution satisfies 

'w + x + y + z= (mod 2), 

as follows by reduction (mod 2). By definition a Descartes quadruple {w,x,y,z) satisfying 
w + x + y + z>Oisa, root quadruple if and only if 

(1) w = —n <0<x<y<z, 

(2) 2{-n + x + y) - z > z. 

Condition (1) implies conversely that w + x + y + z>0, and condition (2) is equivalent to 

—n + x + y — z>0. 

The integer solutions Qviw, x, y,z) =0 with w = —n < are in one-to-one correspondence 
with integer representations of n? by the the ternary quadratic form 

QriX, Y, Z) := XY + XZ + YZ, 

The correspondence between solutions to Q-viw, x, y, z) =0 and solutions to Qt{X, Y, Z) = v? 
is given by 

{X,Y,Z) := {-{w + x + y - z), -{w + X - y + z), -{w - X + y + z)). 

The congruence condition w + x + y + z = Q (mod 2) on integer Descartes quadruples implies 
that (X, y, Z) are all integers. The map is onto, because an explicit inverse map from an 
integer solution (X, Y, Z) is 

(u), X, y, z) := (-n, n + X + Y,n + X + Z,n + Y + Z). 

The primitivity condition gcd{—n, x, y,z) = 1 on Descartes quadruples is equivalent to the 
primitivity condition gcd{X,Y, Z) = 1. The "root quadruple conditions" (1) and (2) above are 
equivalent to the inequalities 

0<X<Y<Z. 

For any integer M, the integer solutions of Qt{X, Y,Z) = M are in one-to-one correspon- 
dence with integer representations of M by the determinant ternary quadratic form 

Qa{A,B,C):=AC-B\ (4.6) 
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A solution [X, Y, Z) gives a solution to Qa.{A, B,C) = M under the substitution 

iA,B,C):=iX + Y,X,X + Z). 

The inverse map is 

{X, Y, Z) ■= {B,A-B,C- B). 

The "root quadruple conditions" above on {X, Y, Z) are easily checked to be equivalent to the 
inequalities 

0<2B <A<C. (4.7) 

In this correspondence. B is necessarily an integer, so 2B is an even integer. The primitivity 
condition gcd{X,Y,Z) = 1 transforms to the primitivity condition gcd{A,B,C) = 1. 

The conditions (|4.7|) give a complete set of equivalence classes of primitive positive definite 
integral binary forms of fixed discriminant D = — 4M under the action of GL(2,Z). To show 
this, note that the conditions for the positive definite binary quadratic form [A, B, C] := 
AT^ + 2BTU + CU"^ of discriminant D = 45^ _ ^AC = -AM to be a reduced form in the 
SL{2,Z) sense are 

< \2B\ <A<C. (4.8) 

It is known that all positive definite integral forms are 5*^(2, Z) -equivalent to a reduced form, 
and that all reduced forms are 5'L(2, Z)-inequivalent, except for [j4,S,A] ~ [^, — S,A] and 
[2A, A, C] ss [2A, —A, C], where ss denotes SL{2, Z)-equivalence. Since GL{2, Z) adds only the 



action of 



1 
-1 



it follows that all are equivalent to a form with -B > 0, and all of these are 
GL(2,Z)- inequivalent. 

Combining these two steps, where one takes M = n'^ in the second step, associates to any 
Descartes quadruple (— n, x, y, z) the definite binary quadratic form 

[A, B, C] := [-n + x,-{-n + x + y-z),-n + y] (4.9) 

of discriminant D := {2B)'^ — AAC = — 4n^. Furthermore this form is reduced in the GL{2, Z)- 
sense ()4.7() if and only if (— n, x, y, z) is a root quadruple. The primitivity condition 
gcd{—n,x,y,z) = 1 transforms to gcd{A,B,C) = 1; conversely, gcd{A,B,C) = 1 implies 
gcd{-n,x,y,z) = 1. D 

Using Dirichlet's class number formula we obtain an explicit formula for N*g^^{—n). 
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Theorem 4.3. Forn > 1 the number N*^^^{—n) of primitive integral root quadruples satisfies 
Kooti-n) = jWi^- ^^^) + 2-W-5.-1, (4.10) 



4 J. J. p 

p\n 



where X-4(f^) = (—1) for n and for even n, uj{n) denotes the number of distinct primes 

dividing n, and (5n = 1 if n = 2(mod 4) and (5n = otherwise. 

Proof. The condition for a positive definite form [A, i?,C] to be 5L(2, Z)-reduced is that 
< 2\B\ < A < C. The form classes under GL{2, Z)-equivalence are derived from the SL{2, Z)- 
equivalence classes using the coset decomposition 



GL{2, 



{I, 



1 
0-1 



}SL{2,Z) 



The action of 



1 

0-1 



maps the set of 5L(2,Z)-reduced forms into itself, hence 

h^{-4n^) = -{h{-4m^) + a{-4m'^)) 



(4.11) 



where a{—4m?) denotes the number of primitive reduced ambiguous forms of discriminant 
— 4?7i^, where we define a reduced ambiguous form to be any S'L(2,Z)-reduced form of the 
shape [A, 0,C], [B,2B,C], and [A,2B,A], with B > 0. (These are exactly the reduced forms 
with S > and [A,B,C] « [A,-B,C], see Mathews (HOI PP- 69-72 ].) A formula for the 
general class number h{d) was given by Dirichlet, and we use the version in Landau [^5]. For 
a discriminant d < 0, Landau 25, Theorem 209], states that 



hid) 



Wd 

2tt 



\d\L{l,Xd) 



in which L{s,Xd) = Yl'^=iXdi'n)n~'^, with Xdin) = (-) a real Dirichlet character (mod \d\), 
and Wd = 2 a d < —4, with w-^ = 4. Thus for d = —4v? with n > 1, 

2n 
/i(-4n2) = — L(l,x_4n2)- 

TT 



Landau j25l Theorems 214] also gives 

^(l,X-4n2) 



p\n 



X-i{p) - 

p ' 



LihX- 
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Combining this with Lambert's formula -L(l,X-4) = j, we obtain, /i(— 4) = 1 and, for n > 1, 
the class number formula 

M-4n^) = fn(l-^)- (4-12) 

p\n 

It remains to determine a(— 4n^). We claim that, for n > 1, 

, ,2a r 2^(")^i if n = 2 (mod 4), ,, ^^, 

ai-An'^) = i „,,/ ^ , . ^ ^' (4.13) 

^ '^ \ 2'^("-' otherwise. ^ ^ 

To prove this we consider for each j > separately the cases n = 2^s, with s odd. There are 

three types of ambiguous reduced forms to consider: 

(1) [A,0,C] withO< A< C; 

(2) [2A, A, C] with 0<2A<C; 

(3) [A,B,A], with < 2B < A. 

The primitivity requirement is that gcd{A, B, C) = 1. For type (1), the discriminant condition 
gives n^ = AC. Since gcd{A, C) = 1 each relatively prime factorization n = pq with gcd{p, q) = 
1 gives a candidate pair {A, C) = {p^ ,q^), and the requirement A < C rules out half of these, 
provided n > 1, for then A = C cannot occur. We conclude that there are exactly 2'^^"''^ 
solutions of this type, when n > 1. The case n = 1 is exceptional, with A = C = 1 as the only 
solution. 

For the remaining count we must group types (2) and (3) together. For type (2) we have 
-4n^ = 4A'^ - SAC which gives n^ = A{2C - A). The condition gcd{A,C) = 1 implies 
gcd{A, 2C - ^) = 1 or 2. For type (3) we have -4n^ = 45^ - 4A'^ which gives n'^ = {A + 
B){A - B). The condition gcd{A, B) = I implies gcd{A + B,A - B) = 1 or 2. 

We first suppose j = l- This is the easiest case because there are no solutions of types (2) 
and (3). Indeed, for type (2) we have A{2C — A) is even, so A must be even, and 2C — A 
is even. But gcd{A,C) = 1 means C is odd, so A and 2C — A are not congruent (mod 4). 
Thus one of them is divisible by 4, hence 8 divides A{2C — A) = n? which contradicts v? being 
divisible by 4 but not 8. For type (3) we have {A + B){A — B) is even, so both factors are even. 
At least one of A, B is odd since gcd{A,B) = 1, so the two factors are incongruent (mod 4), 
hence one of them is divisible by 4, so again we deduce that 8 divides n^, a contradiction. 

We next suppose j = 0. Then n is odd and the primitivity requirement for type (2) forms 
becomes gcd{A, 2C—A) = 1, and for type (3) forms becomes gcd{A+ B , A — B) = 1. Consider a 
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factorization r? = p^q^ with gcd(p, q) = 1. For type (2), setting this equal to the factorization 
A{2C - A) yields A = p'^,C = l/2(p2 + q'^). The condition 2A < C becomes 

2 ^ -l- 2 

For type (3), setting n = q^p^ equal to the factorization {A + B){A — B) yields A = l/2{q^ + 
p^), B = \/2{q^ — P^)- The condition B >Q gives p^ < q^ and the condition 2B < A gives 
3p^ > g2; thus g^/3 < p^ < q^, and we have g^/3 7^ p'^ by the relative primality condition. 
These two cases both require p'^ < q^, and they have disjoint ranges for p^ and cover the whole 
range < p^ < q^- Thus exactly half the factorizations lead to one solution and the other half 
to no solution, and we obtain exactly 2'^'"' ""^ solutions of types (2) and (3). 

Finally we suppose j > 2. In these cases the primitivity requirement for type (2) forms 
becomes gcd{A, 2C — A) = 2, and for type (3) forms becomes gcd{A + B,A — B) = 2. We 
have n^ = 2'^^p'^q^ and there are several splittings of factors to consider, for example A = 
2^^- V, 2C - A = 2q^ OT A = 2q^ ,2C - A = 2^^-^p^ in type (2) andA + B = 2^^-^q'^, A-B = 
2p^ or vice versa in case (3). In all cases, a careful counting of factorizations finds exactly half 
the possible factorizations lead to solutions of types (2) and (3), giving 2^^'^'^ solutions of 
types (2) and (3). We omit the details. 

Adding up these cases gives 2'^^^> solutions in total when j = or j > 2 and 2^^"'>^ 
solutions if j = 1. This proves the claim. 

Combining (|4.11|) and (|4.12|) with the claim yields the desired result. □ 

The exact formula for N*^g^{—n) for prime n = p gives N*^^^{—p) = 2i- if p = l(mod 4) 
and N*g^f{—p) = 2l- if p = 3(mod 4), as in Table 2. The exact formula also yields the 
following asymptotic upper and lower bounds for N*^^^{—n). 

Theorem 4.4. There are positive constants Ci and C2 such that, for all n>2>, 



n 
log log n 



CiTTTTTrr < Kooti-n) <C2n log log n. (4.14) 



Proof. We use the exact formula H4.1U|) of Theorem 14.31 The second term ia(— 4n^) grows 
like 0{n'^) and is negligible compared to the first term. To estimate the first term we use 



27 



inequalities of the form 



and 



-'■-'■ p log log n 

p\n 



TT(1 + -) < csloglogn. (4.16) 

P 

p\n 

To establish ()4.16|) . we note that its left side is upper bounded by the product taken over the 
first k distinct primes, with k chosen minimal so that the product exceeds n. Thus 

log(JJ(l + -)) < Yl log(l + -)< Yl -<logloglogn 

p\n p<log n/ log log ra p<log n/ log log n 

for large enough n, using 20, Theorem 427], and exponentiating this yields (|4.16|) . The in- 
equality H4.15() follows from this because Yini^ j) converges. □ 

The bounds of Theorem 14.41 are sharp up to a multiplicative constant. We can take n 
to be a product of all small primes p = 3 (mod 4) to achieve a sequence of values with 
A^*QQj(— n) >> n log log n, and to be a product of all small primes p = 1 (mod 4) to achieve a 
sequence of values with N;^^^{-n) « i^p^. 

Remark. Theorem 14 . 21 gives a one-to-one correspondence between Descartes quadruples con- 
taining —n and integral binary quadratic forms of discriminant —An? having nonnegative mid- 
dle coefficient. These Descartes configurations are exactly the ones that contain the outer 
circle of the packing as one of their circles. Under this correspondence the reduction algorithm 
of §3 for Descartes quadruples matches the Gaussian reduction algorithm for positive definite 
binary quadratic forms. 

A related correspondence holds for the exceptional integral Apollonian packing (0,0,1,1) 
having n = 0. There is a unique primitive reduced positive definite form of determinant 0, 
namely Q = [0,0,1] = C/^, which is an ambiguous form. In this case the special Descartes 
quadruples are those that contain one of the two infinite edges of the packing as one of their 
circles. They are essentially Descartes configurations of "Ford circles", and on these config- 
urations the reduction algorithm corresponds to an additive variant of the continued fraction 
algorithm. For futher information see Ford J2j or Rademacher j38l pp. 41-46], |39l pp. 
264-267]. 
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5. Integers Represented by a Packing: Asymptotics 

In this section we study the ensemble of integer curvatures that occur in an integer Apohonian 
circle packing, where integers are counted with the multiplicity that they occur in the packing. 
Their asymptotics are known to be related to the Hausdorff dimension of the residual set of 
the packing, as follows from work of Boyd described in Theorem 15.21 below. At the end of 
the section we begin the study of the set of integer curvatures that occur, counted without 
multiplicity. 

The residual set of a disk packing V (not necessarily an Apollonian packing) is the set 
remaining after all the (open) disks in the packing are removed, including any disks with 
"center at infinity." For a general disk packing V, we denote the Hausdorff dimension of the 
residual set by a{V) and call it the residual set dimension of the packing. The definition 
of Hausdorff dimension can be found in Falconer ^1], who also studies the residual sets of 
Apollonian packings in J14l pp. 125-131.]. 

The residual sets of Apollonian packings all have the same Hausdorff dimension, which 
we denote by a. This is a consequence of the equivalence of such residual sets under Mobius 
transformations (see jl71 Sect. 2]), using also the fact that the Hausdorff dimension strictly 
exceeds one, as follows from results described below. 

The exponent or packing constant e(V) of a bounded circle packing V (not necessarily an 
Apollonian packing) is defined to be 

e{V) := sup{e : ^ r{Cy = oo} = inf{e : ^ r{Cy < cx)}, 
cev cev 

in which r{C) denotes the radius of the circle C. This number has been extensively studied 

in the literature, beginning in 1966 with the work of Melzak |33| Theorem 3], who showed 

that in any circle packing that covers all but a set of measure zero one has X^ceP '"(^) ~ ^^■ 

He constructed a circle packing with e('P) = 2 and showed for Apollonian packings that e(P) 

lies strictly between 1.035 and 1.99971. He conjectured that the minimal value of e{V) is 

attained by an Apollonian circle packing. In 1967 J. Wilker ^;51j showed that all osculatory 

circle packings V, which include all Apollonian circle packings, have the same exponent e{V), 

which we call the osculatory packing exponent e. He also showed that e > 1.059. Later Boyd 

[H]i 1^1 |Zj improved this to 1.300 < e < 1.314. Recent non-rigorous computations of Thomas 
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and Dhar 1^ estimate the Apollonian packing exponent to be 1.30568673 with a possible error 
of 1 in the last digit. 

The relation between the packing exponent and the residual set dimension of Apollonian 
packings was resolved by an elegant result of D. Boyd jQ. 

Theorem 5.1. (Boyd) The exponent e of any bounded Apollonian circle packing is equal to 
the Hausdorff dimension a of the residual set of any Apollonian circle packing. 

The inequality e > a follows from a 1966 result of Larman ^7\, and in 1973 Boyd proved 
the matching upper bound a > e. A simpler proof of the upper bound was later given by C. 
Tricot gnj. 

Given a bounded circle packing V we define the circle- counting function N-p{T) to count 
the number of circles in the packing whose radius of curvature is no larger than T, i.e., whose 
radius is at least y. Boyd i7_ proved the following improvement of the result above. 

Theorem 5.2. (Boyd) For a bounded Apollonian circle packing V , the circle- counting func- 
tion N'p{T) satisfies 

T^oo logT ' ^ ^ 

where a is the Hausdorff dimension of the residual set. That is, N'p{T) = 7""+o(i) as T ^ oo. 

. Theorem 13.31 showed that the curvatures of all circles in the packing, excluding the root 
quadruple, can be enumerated by the elements of the Apollonian group A. ^From this one 
can derive a relation between the number of elements of A having height below a given bound 
T and the Hausdorff dimension a. We measure the height of an element M £ A using the 
Frobenius norm 

||M||f := (ir[M^M])i/2 ^ (^^2y/2_ ^5^2) 

id 

Theorem 5.3. The number of elements Nt{A) of height at most T in the Apollonian group 
A satisfies 

A^t(^) = T°+°«, (5.3) 

as T —> oo, where a is the Hausdorff dimension of the residual set of any Apollonian packing. 
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In order to prove this result, we establish two preliminary lemmas. 

Lemma 5.4. Let M = Si^ ■ ■ ■ Si^Si-^ G A, the Apollonian group, and suppose that ij ^ ij^i 
for I < j < m — 1, and m> 2. In each row k of M , 

(i) Mki<0 ifl = h, 

(a) Mkj > \Mki\ for I = h and j ^l. 

Proof. The lemma follows by induction on m. It is true for ttt, = 1, since each matrix Si has 
i column negative (or zero). 

Suppose (i)-(ii) hold for M' = Si^ ■ ■ ■ Si^. Suppose, for convenience, that ii = 1. Then 



M = M'Si- 



-M[^ 2M[^ + M{2 2M[^ + M[^ 2M[^ + M[^ 

-M'^^ 2M'^^ + M32 2M;^i + M;^3 2M31 + M34 

-Mil 2Mii + Mi2 2Mii + M^3 2Mii + M'^^ j 



Since ^2 / ^1 = 1 all M[^ > by (ii) of the induction hypothesis, so Mn = M[^ < gives (i). 
Next, note that 

Mkj = 2M',, + M',^ > 2M',, -\M',i\> M',, = \Mki\ 

since ML > |ML| and ML- > — |M^J in all cases by (ii). This completes the induction step in 
this case. The arguments when ii = 2, 3, or 4 are similar. □ 

Lemma 5.5. Let v = (a, b, c, d) be an integer root quadruple with a < 0. Then there are 
positive constants cq = cq(v) and ci = ci(v) depending on v such that 



co||M||ir < ||Mv||oo < ci||M||i7, for all Me A. 



(5.4) 



Proof. For the upper bound, we have 



||Mv||oo < 2|Mv| < 2||M||f|v| , 



(5.5) 



so we may take ci = 2|v|. 

For the lower bound, we first show that if M = Si^ 
have 

||Mv||oo>^||M||p 



■ Si^Si-^ with ij / ij^i and ii = 1, we 



(5.6) 
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The vector v has sign pattern (—,+,+,+) and Lemma |5.4I shows that M has first column 
nonpositive elements and other columns nonnegative. Thus all terms in the product Mv are 
nonnegative, and hence 

4 4 

(Mv)i > Y^ \Mij\\-Vj\ >Y\Mij\, 
because a < implies min(|o|, |6|, |c|, \d\) > 1. Thus 

l|Afv||oo > -y^lMjJ > -\\M\\f . 

M MOO _ ^ ^_^ I iji _ 2 II ll-f* 

It remains to deal with the cases where ii = 2, 3 or 4. By Theorem 13 .31 the value ||Mv||oo 
gives the curvature of a particular circle in the packing, and this circle lies in one of the four 
lunes pictured in Figure 3 according to the value of im ■ 




Figure 4: Four lunes of Descartes quadruple. 

The bound ()5.6() applies to all circles in the central lune corresponding to im = 1- For 
the remaining cases, we use the fact that there exists a Mobius transformation : C ^ C 
with (j) G Aut{V), which fixes the Descartes configuration corresponding to v but cyclically 
permutes the four circles a ^>- b ^>- c ^ d. In particular (f) also cyclically permutes the four 
lunes zi = 1 ^ zi = 2 — > ii = 3 — > ii = 4. Now (p maps the center of circle d to the center of 
circle a, which is the point at infinity, and maps the point at infinity to the center of circle b. 
It follows that the stretching factor of the map <j) inside the four lunes is bounded above and 
below by positive absolute constants C2 and c^ . Since 4> maps the lune ii = 4 to ii = 1 we 
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conclude for cases where ii = 4 that 

||Mv||oo>;^||M||^ . (5.7) 

ZC2 

Applying the same argument to 4>^ and (f)^ gives the similar bound for the cases ii = 3 and 
ii = 2. We conclude that the lower bound in (|5.4() holds with cq = i^- □ 

Proof of Theorem 15.31 Pick a fixed quadruple having a < 0, say v = (—1,2, 2, 3), and let 

Pv he the associated Apollonian packing. By Theorem l3.3l each M £ A corresponds to a circle 
of curvature ||Mv||oo in V^, and all circles are so labelled except the four circles in v. Lemma 
15.51 shows that each ||M||i? < T produces a circle of curvature at most ciT. Now Theorem 15.21 
asserts there are at most T°'+°(^) such circles, hence Nt{A) < T°'+°(^) , Conversely Lemma 
15.51 implies that each circle of curvature ||Mv||oo < T comes from a matrix M £ A with 
II-^IIf ^ —T. Since there are at least 7"°+°(i) such circles, we obtain Nt{A) > r"^°'^\ as 
desired. □ 

Can the estimate of Theorem 15.31 be sharpened to obtain an asymptotic formula? A. 
Gamburd has pointed out to us that the method of Lax and Phillips 28 might prove useful 
in studying this question. 

We now turn to a different question: How many different integers occur, counted without 
multiplicity, in a given integral Apollonian circle packing Vv? This seems to be a difficult 
problem. It is easy to prove that at least cT^''^ of all integers less than T occur in a given 
packing. This comes from considering the largest elements of the vectors {(5'iS'2)-'v : j = 
1,2, ...}, where v is a root quadruple, which are curvatures in the packing, by Theorem 13.31 
above. These values grow like j^ (see the example (1) in §7). Concerning the true answer to 
the question above, we propose the following conjecture. 

Positive Density Conjecture. Each integral Apollonian packing represents a positive frac- 
tion of all integers. 

Theorem 15.21 shows that the average number of representations of an integer n grows like 
n"~^, which goes rapidly to infinity as n ^ oo. Therefore one might guess that all sufficiently 
large integers are represented. However in the next section we will show there are always some 
congruence restrictions on which integers occur. There we formulate a stronger version of this 
conjecture and present numerical evidence concerning it. 
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6. Integers Represented by a Packing: Congruence Conditions 

In this section we study congruence restrictions on tlie set of integer curvatures whicli occur 
in a primitive integral Apollonian packing. 

We first show that there are always congruence restrictions (mod 12). 

Theorem 6.1. In any primitive integral Apollonian packing, the Descartes quadruples (mod 12) 
all fall in exactly one of four possible orbits. The first orbit Y (mod 12) consists of all permu- 
tations of 



{(0,0,1,1), (0,1, 1,4), (0,1, 4, 9), (1,4, 4, 9), (4, 4, 9, 9)} (mod 12). 
The other three orbits are (3, 3, 3, 3) - Y, (6, 6, 6, 6) + F and (9, 9, 9, 9) - F (mod 12). 



(6.1) 



Remark. Each orbit contains only 4 different residue classes (mod 12), hence 8 residue classes 
(mod 12) are excluded as curvature values. 

Proof. A straightforward computation, using the action of the Apollonian group (mod 
12), shows that (up to ordering) the set of all quadruples without common factors of 2 or 3 
(mod 12) consists of the list below, which are grouped into eight orbits under the action of the 
Apollonian group (mod 12). 



(1) 


Y 


(0,0,1,1) 


(0,1,1,4) 


(0,1,4,9) 


(1,4,4,9) 


(4,4,9,9) 


(2) 


(3,3,3,3) -y = 


(6,6,11,11) 


(2,6,11,11) 


(2,3,6,11) 


(2,2,3,11) 


(2,2,3,3) 


(3) 


(6,6,6,6)+y = 


(3,3,10,10) 


(3,6,7,10) 


(3,7,10,10) 


(6,7,7,10) 


(6,6,7,7) 


(4) 


(9,9,9,9) -y = 


(0,0,5,5) 


(0,5,5,8) 


(0,5,8,9) 


(5,8,8,9) 


(8,8,9,9) 


(5) 


-y 


(0,0,11,11) 


(0,8,11,11) 


(0,3,8,11) 


(3,8,8,11) 


(3,3,8,8) 


(6) 


(3,3,3,3)+y = 


(0,0,7,7) 


(0,4,7,7) 


(0,3,4,7) 


(3,4,4,7) 


(3,3,4,4) 


(7) 


(6,6,6,6)-y = 


(2,2,9,9) 


(2,2,5,9) 


(2,5,6,9) 


(2,5,5,6) 


(5,5,6,6) 


(8) 


(9,9,9,9)+y = 


(9,9,10,10) 


(1,9,10,10) 


(1,6,9,10) 


(1,1,6,10) 


(1,1,6,6) 



To check the orbit structure is as given, note that the action of the four generators of the 
Apollonian group on the five elements of the orbit Y is summarized in the following transition 
matrix: 





(^ 


2 








\ 


1 

4 


1 


1 


2 











1 


2 


1 











2 


1 


1 




^0 








2 


2/ 
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We may view this matrix as the transition matrix of a Markov chain (after rescaUng each 
row to be stochastic), and find that the action is transitive and the stationary distribution is 
(to'5'5'5'To)' "^^^ other seven orbits have the same transition matrix and the same stationary 
distribution as Y. 

There exist integral solutions to the Descartes equation in all of the congruence classes 
(mod 12) in the list above. However we recall that a Descartes quadruple (a, b, c, d) coming 
from an Apollonian packing satisfies the extra condition 



a + b + c + d>0. 



(6.2) 



In the rest of the proof we show that this extra condition excludes half of the orbits above, 
namely orbits (5)- (8). 

As a preliminary, we observe that any integer solution (a, b, c, d) to the Descartes equation 
Hl.l|) yields a unique integer solution to the equation 



2 I ^2 _ ;2 



4m^ + 4a^ + n^ = / 



(6.3) 



and vice-versa. Here the solution to 1)6. 3|) is given by 



a 

n 

I 

m 



1 

1 

2 1 

1 _i 

'2 2 




-1 

1 

_1 1 

2 2 J 



a 
b 
c 
d 



a 

b-c 

2a + b + c 

^{d — a — b — c) 



(6.4) 



In the reverse direction, an integer solution to (|6.3|) gives one to the Descartes equation via 



a 
b 
c 
d 



10 



-1 



1 1 

2 2 
1 1 





- 2 2 

-10 12 



a 

n 

I 

m 



a 

\{l -2a + n) 

\{l-2a-n) 

2m + I — a 



(6.5) 



Solutions to the Descartes equation satisfy a congruence (mod 2) which guarantee that the 
maps above take integral solutions to integral solutions, in both directions. Now (|6.3|) gives 



f > Aa^ + m^ > 2(\a\ + \m\f > (\a\ + |m|)^ 



(6.6) 



and equality holds if and only if £ = a = m = 0. In particular, if ^ > 0, then (|6.6j) gives 



> \a\ + \m\. 



(6.7) 
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We assert that any integer solution (a, 6, c, d) to the Descartes equation has 

a + b + c + d>0 if and only if / > 0. (6.8) 

To prove this, note that ifa-\-b + c + d>0 then by Lemma l3. II fi) we have 

e = 2a + b + c = {a + b) + {a + c) > 0. 

Equality can hold here only if a = 6 = c = 0, which implies d = 0, which contradicts the 
assumption a + b + c + d > . Conversely, if £ > 0, then, using (|6.7|1 . 

a + b + c + d = 2i + 2m-2a>2{£-\a\- \m\) > 0, 

so (|6.8() is proved. 

Claim : No primitive integer Descartes quadruples with a + b + c + d > occur in the 
orbits (5)-(8). 

We prove the claim for orbit (8); the arguments to rule out orbits (5), (6), (7) are simi- 
lar. We argue by contradiction. Suppose there were such a solution in orbit (8). Since the 
Apollonian group acts transitively on the orbit, and preserves the condition a + b + c + d > 0, 
there would be such a quadruple (a, b, c, d) = (1, 1, 6, 6)(mod 12). In this case l = 2a + b + c = 
9(mod 12) and 

m=-(6 -6 + 1 + 1) = l(mod 6), 

which gives m? = l(mod 12). Now (|6.3jl gives 

(/ + 2m){l - 2m) = l^ - Arn^ = Aa^ + n^ > 0. (6.9) 

Since a + b + c + d > we have I > hy ()6.8|) . Then in the equation above at least one of the 
factors on the left side must be positive, hence they both are. Consider / + 2m > 0. We have 

^ + 2m = 9±2(mod 12) = 3(mod4). (6.10) 

Consider any prime p = 3(mod 4) dividing / + 2m,. Then it divides 4a^ + n^, which it must 
divide to an even power, say p^"^, with a = n = O(mod p'^). If p also divides / — 2m,, then 
it would divide both / and m, and then (|6.5|) would imply that it divides gcd(a, b, c, d), 
which contradicts the primitivity assumption gcd(a, b, c, d) = 1. Therefore p does not divide 
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/ — 2m, and we conclude from ()6.9() that p^^Wl + 2m. It follows that all primes p = 3(mod 4) 
that divide / + 2m do so to an even power, hence we must have have / + 2m = l(mod 4), a 
contradiction. This rules out orbit (8), which proves the claim in this case. 

Theorem 16.11 follows from the claim. □ 

At the end of this section we present numerical evidence that suggests that these con- 
gruences (mod 12) are the only congruence restrictions for the integer packing (—1,2,2,3). 
However there are stronger modular restrictions (mod 24) that apply to other integer pack- 
ings. For example, in the packing (0,0,1,1) (Fig. 2), any curvature which occurs must be 
congruent to 0,1,4,9,12 or 16 (mod 24) (these are the quadratic residues modulo 24). Thus 
only 6 classes (mod 24) can occur rather than the 8 classes allowed by Theorem 16. 11 

It seems likely that the full set of congruence restrictions possible (mod m) is attained for 
m a small fixed power 2^3'', perhaps even m = 24. We are a long way from proving this. As 
evidence in its favor, we prove the following result, which shows that all residue classes modulo 
m do occur for any m relatively prime to 30. 

Theorem 6.2. Let V he a primitive integral Apollonian circle packing. For any integer m 
with gcd{m, 30) = 1, every residue class modulo m occurs as the value of some circle curvature 
in the packing V. 

Proof. Observe that the s-term product W{s) = . . . S2S1S2S1 is 

/ -s s + 1 s(s + l) s{s + l) \ 

-(s-1) s s(s-l) s{s-l) 

10 

V 1 / 

(where the top two rows are interchanged if s is even). Of course, the two non-trivial rows can 
be placed anywhere by choosing the two matrices from the set Si, 52, 5*3, 84^ appropriately. 
If (a, b, c, d)'^ is a quadruple in V then the product 

W{s){a,b,c,df = 
{-sa + {s + 1)6 + s{s + l)c + s{s + l)d, -{s - l)a + sb + s{s - l)c + s{s - l)d, c, d^ (6.11) 

is also in V as well. Let J denote the set of all rows (q, /3, 7, 5) which can occur in a product 
of matrices taken from ^i, 5*2, 5*3, ^4. Thus, if (a, /J, 7, 5) ^ J then so are: 
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(-a, 2a + /3, 2a + 7, 2a + 5), (2/3 + a, -/?, 2/3 + 7, 2/3 + (5), (27 + a, 27 + /3, -7, 27 + 5), 

and (25 + a, 2<5 + P, 25 + 7, -6). Therefore, 

(-a, 2a + /3, 2a + 7, 2a + 5), (3a + 2/3, -2a - /3, 6a + 2/3 + 7, 6a + 2/3 + (5), 

(—3a — 2/3, 4a + 3/3, 12a + 6/3 + 7, 12a + 6/3 + 5), . . ., and in general, 

{-ra - (r - l)/3, (r + l)a + r/3, r(r + l)a + r{r - l)/3 + 7, r(r + l)a + r{r - l)/3 + 6) (6.12) 

are all in J for all r (as well as all permutations of these). Now substitute (a,/3,7,(^) = 
{s{s + 1), s{s + 1), — s, s + 1) G J" into (|6.12j) . This shows that the row vector 

p ■= (_(2r - l)s(s + 1), (2r + l)s(s + 1), 2r2s(s + 1) - s, 2r^s{s + l)+s + l)£j. (6.13) 

The sum of the last two coordinates of p is 

Ar'^s{s + 1) + 1 = r2((2s + 1)^ - 1) + 1 = r^{x'^ _ 1) + 1 = ^^ - r^ + 1 (6.14) 

where u = rx and x = 2s + 1. The g.c.d. of these two summands must divide their difference, 
which is 2s + 1. It is well known (and easy to show) that for any prime power p'^ with p > 5, in 
at least one of the pairs {1, 2}, {4, 5} and {9, 10} are both nonzero quadratic residues modulo 
p^. For each p \ m, let {ap,ap + 1} denote such a pair. Define Up and r^ so that 

Up = Up, Tp = Op + 1 (mod p^) (6.15) 

where p"^ is the largest power of p dividing m. Since gcd{rp,p) = 1 then we can define 
Xp = UpVp^ (mod p"^). We can guarantee that Xp is odd by adding a multiple of p^i if necessary. 
Hence, for these choices, the expression in ()6.14|) is modulo p'^, i.e., 

r^(x^-l) + l = (modp'"). (6.16) 

Of course, we can use the values Vp + kp^ and Xp + Ip"^ in place of r^ and Xp in 1)6. 16() for any 
k and /. Note that gcd(xp,p) = 1. Letting p range over all prime divisors of m, then by the 
Chinese Remainder Theorem, there exist X (odd) and R such that 

R\X^ - 1) + 1 = (mod p"") (6.17) 

for all p'^l m. Thus, 

R^iX'^ -1) + 1 = (mod m), gcd{X, m) = 1. (6.18) 
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Hence, by ()6.13|) and H6.14() we can find a row modulo min J' of the form (C, D, A, — ^)(mod m) 
where it easy to check that gcd{A,m) = 1. 

We can now apply the transformation preceding (|6.12j) to (C, D, A, —A) to get the following 
rows modulo m in J: 

{C, D, A, -A) (modm) 

{2A + C, 2A + D, ~A,A) (modm) 

{AA + C, AA + D, A, -A) (modm) 

(6^ + C, 6A + D, -A, A) (mod m) 

and more generally 

{itA + C, UA + D, A, -A) (mod m) £ J for ah t > 0. (6.19) 

Suppose for the moment (and we will prove this shortly) that we can find (a, b, c, d) S V with 
gcd(a + 6, m) = 1. Taking the inner product of the row in H6.19() with (a, 6, c, d)'^ , we get the 
curvature value 

{AtA + C,4:tA + D,A,-A)-{a,h,c,dY (modm) (6.20) 

= AA{a + h)t + Ca + Db + Ac- Ad (modm). (6.21) 

Since gcd(4A(a + b),m) = 1 then these values range over a complete residue system modulo 

m as t runs over all positive integers. The proof will be complete now if we can establish the 

following result. 

Claim. If "P is a primitive packing then for any odd m > 1, there exists (a, b, c, d) G V with 

gcd(a + 6, m) = 1. 

Proof of Claim. First recall that for any (a,b,c,d) £ V, we have gcd(a,6, c) = 1. We have 

also seen by ()6.11() . if (a, b, c,d) £ V then for any r > 1, 

(A(r),i?(r),C(r),D(r)): = 
(-ra+(r + l)6 + r(r + l)(c + (i),-(r- l)a + rb + r{r - l){c + d),c,d) £ V (6.22) 

as well. Define q{r) to be the sum of the first two components of this vector: 

q{r) ■= A{r) + B{r) = 2(c + dy"^ - 2{a -b)r + a + b. 
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Let p denote a fixed odd prime. We show that 

q{r) ^ (mod p) for some r > 1. (6.23) 

Suppose to the contrary that q{r) = (mod p) for all r. Thus, 

g(0) = a + 6 = (mod p) 

q{l) = 2(c + d) -2(a-5) + (a + 6) =2(c + d) -a + 36 = (modp) 

q{2) = 8(c + d)- 4(a - 6) + (a + 6) ee 8(c + d) - 3a + 56 = (mod p) 

which implies a = 6 = c + (i = (mod p). However, since 

a = b + c + d^ 2y^6(c + d) + c d then cd = (mod p), i.e., c = or d = (mod p). 

This would imply that V is not primitive, a contradiction which establishes ()6.23|) . 

To finish proving the claim, for each p\m let rp satisfy q{rp) ^ (mod p). Then we have 

q{rp + kp) = q{rp) ^ (mod p) 

for all k > 0. By the Chinese Remainder Theorem one can find R and S such that q{R + kS) ^ 
(mod p) for all p | m and all k. In particular, gcd{q{R),m) = gcd{A{R) + B{R),m) = 1, and 
the Claim is proved. □ 

Which integers occur as curvatures, when the congruence conditions are taken into account? 
We consider numerical data for two cases. The first case is the packing with root quadruple 
(—1,2,2,3), where Theorem 16.11 permits only values 2,3,6, or 11 (mod 12). Not all such 
integers appear in the Apollonian packing (—1,2,2,3), for example in the class 6 (mod 12) 
the value 78 is missed. In Table 3 we present the missing values in these residue classes for 
the first million integers. Only 61 integers congruent to 2, 3, or 6 do not occur in the packing 
(—1,2,2,3), the largest being 97287 (see Table 3), and no integers 11 (mod 12) are missed. 
This data suggests that there are finitely many missing values in total, with 97287 being the 
largest one. 

Our second example is the packing with root quadruple (0,0,1,1). As mentioned above, 
there are congruence conditions (mod 24) in this case. Table 4 presents numerical data on 
exceptional values for the allowed congruence classes (mod 24) up to T = 10^. There is a 
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n = 3 (] 


mod 12) 










159 


207 


243 


435 


603 


711 


1923 


2175 


2319 


3711 


4167 


4959 


4995 


5283 


6015 


6879 


7863 


10095 


10923 


11295 


12063 


16311 


16515 


18051 


19815 


21135 


23175 


28323 


41655 


48075 


68055 


97287 



















n = 6 (mod 12) 



78 


246 


342 


834 


1422 


2010 


2022 


2454 


2718 


2766 


3150 


3402 


3510 


3774 


4854 


6018 


6666 


7470 


10638 


12534 


13154 


13206 


20406 


24270 


32670 


42186 


45258 


55878 







n = 2 (mod 12) 



13154 



Table 3: Missing integers in the packing (—1,2,2,3) up to 10^ 

much larger set of exceptional values, and it appears more equivocal whether the full list of 
exceptional values is finite. However we think it is. 

The numerical examples above support the idea that for any fixed integer Apollonian 
packing and for sufficiently large integers a finite list of congruence conditions will be the 
only obstruction to existence. We therefore propose the following strengthening of the Density 
Conjecture. 

Strong Density Conjecture. In any primitive integral Apollonian packing, all sufficiently 
large integers occur, provided they are not excluded by congruence conditions. 

In further support of the Strong Density Conjecture, we note an analogy to a number- 
theoretic conjecture of Zaremba |S3], who conjectured that there exists an absolute constant 
b (possibly b = 5) such that each sufficiently large positive integer can be represented by some 
continuant with digits bounded above by b. In other words, given any integer m > 1, there 
exists an integer a < m {a relatively prime to m) such that the simple continued fraction 
[0,ci,--- ,Cr] = a/m has partial denominators Cj < b. Fix h, and let M be the set of all 
pairs (o, m) with the above property. There is a linear recurrence for the pairs (o, m) which is 
similar to that of the Descartes quadruples, since if the terms in the continued fraction of a/m 
are bounded by 6, then so are those for the fractions l/{i + a/m), i = 1,2, ... ,b. Zaremba's 
conjecture is saying that all the integers m will appear in some pair of M. This conjecture is 
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24) 








48 


120 


360 


528 


552 


720 


888 


912 


1080 


1176 


1272 


1392 


1560 


1704 


1848 


1968 


2184 


2208 


2736 


2880 


3240 


3408 


3552 


4080 


4392 


4464 


4584 


4680 


4896 


5040 


5088 


5760 


6192 


6888 


7272 


8280 


8880 


9792 


10680 


10920 


10944 


11760 


11928 


13152 


14160 


14328 


16008 


17160 


17232 


17520 


18000 


19320 


20712 


23160 


25896 


26472 


26760 


27552 


27600 


27768 


29424 


29688 


30288 


31440 


34440 


34488 


35232 


36408 


36648 


36816 


37968 


38928 


39168 


43056 


43392 


45240 


46056 


50448 


52800 


58728 


59400 


66120 


74976 


80280 


82200 


87192 


93216 


96912 


96960 


107016 


108240 


117480 


121680 


133392 


137280 


138360 


165360 


201480 


399000 


424560 


496080 























n = 


= 12 (mod 


24) 








132 


252 


300 


468 


636 


780 


1140 


1476 


1572 


1980 


2100 


2148 


2628 


2820 


2868 


3012 


3492 


3828 


3900 


4212 


4692 


5028 


5148 


5340 


5796 


6516 


6684 


6900 


7380 


7908 


8772 


10020 


10212 


10260 


10380 


10548 


11268 


11868 


12876 


13572 


14100 


14244 


14724 


14916 


15300 


15588 


19260 


19620 


20940 


21732 


22908 


23652 


24252 


24804 


25140 


25812 


26100 


26124 


27660 


28860 


29532 


30540 


31092 


31932 


36564 


37908 


38772 


39780 


41460 


41964 


44988 


46980 


52260 


52788 


61596 


67308 


69324 


69420 


75900 


76908 


79740 


88140 


101940 


120300 


135252 


185580 


188748 


220308 


228780 


234660 


354540 


422820 


472548 


926820 


1199820 















n = 1,4 or 9 (mod 24) 








241 


340 748 2980 5452 11380 45652 


16617 21825 






n = 16 (mod 24) 




208 328 712 1168 2488 3400 5200 13600 15088 116896 



Table 4: Missing integers in the packing (0,0, 1, 1), up to 10^ 
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currently still open. But as in the Apollonian packing, consideration of the Hausdorff dimension 
of the set E^ = {a/m : (a, m) G M} is suggestive. Namely, let Sb{m) be the number of a's 
such that {a,m) G M. If Sb{m) ~ m^ , then ^m'^JTi"^ converges iff x > /? + 1. Since the 
abscissa of convergence of the series Yl Sh{ra)m~^ is equal to twice the Hausdorff dimension 
7 of Eb (see T. Cusick [10]), then /3 = 27 - 1 ss .0624 > 0. Thus the "expected" number of 
appearances of m in the pairs of M is m^ ^ 1. 

7. The Growth of Descartes Quadruples in a Packing 

The circles in an integral Apollonian circle packing, starting from the root quadruple, are 
enumerated by the elements of the Apollonian group. The graph of this group is a rooted 
infinite tree with four edges meeting each vertex, with each vertex labelled by a nontrivial 
word in the generators of the Apollonian group. (Such a word satisfies the condition that any 
two adjacent generators in the word are unequal.) Starting from the root node, there are 4 
nodes at depth 1, and at each subsequent level there are three choices of generators at each 
node, so there are 4 x 3"~^ words of length n labelling depth n circles. How are the curvatures 
of the circles at depth n distributed? We consider the maximum value, the minimum value, and 
the median value. In the process we also determine the joint spectral radius of the generators 
of the Apollonian group. 

We begin with the maximum value. We define for n = 4m + i with < i < 3, the reduced 
word T„ of length n given by 

r„ := Ti{SiS3S2Sir, (7.1) 

with Ti = I,Si, S2S1, S3S2S1 for < i < 3, respectively. 

Theorem 7.1. Let v = {a,b,c,d) be any root quadruple with a < b < c < d and a < 0, 

a + b + c + d > 0. Then for any reduced word W of length n in the generators {Si, S2, S3, S4} 
of the Apollonian group, 

||Vrv||oo<||r„v||oo . (7.2) 

Proof. Write W = Si^Si^_-^ ■ ■ ■ Si-^ and set w^") = T^v and v^") = r„v. Write the elements of 
w^"^ and v'"^ in increasing order as 

^(") < w^^^ < w^^^ < w^^^ and ?;i") < v^^^ < v^^^ < v^^^ . 
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The idea of the proof is that T„ always inverts with respect to the circle of smallest curvature, 
and in fact produces the largest curvature vector in a strong lexicographic sense. More precisely, 
we prove by induction on n > 1 that 

Jn) ^ ^(n) for 1 < i < 4 (7.3) 

and 

(n) (n) ^ (n) (n) /rj ^\ 

wl — wl < vl — vl . (7.4) 

For the base case n = 1, we have v'--'^^ = (a', b, c, d) where a' = 2(6 + c + d) — a = | jS'ivl |oo- If 
b' = 2(a + c + d) — b = ||S'2v||oo and c' = HS'svUoo, d' = ||54v||oo then a < b < c < d gives 
d' < c' <b' < a' , and (|7.3() holds for n = 1, since d' > d because v is a root quadruple. 

For the induction step, a reduced word has i„ 7^ in-i- The forbidden move Si^_-^ is the one 
that replaces wl^ with 2{w\^ +^2 "'"'"^3 ) ~ ^4 • Now the induction hypothesis 
gives 

(n) , (n-1) , (n-1) (n) 

W\ < 1(72 < ^2 = V\ 

(n) , (n-1) , (n-l) (n) 

W2 < "W^a ^ ^3 = ^2 

(n) , (n-1) , (n-1) (n) 

fi'a < li'4 < ^"4 = ^^3 



and 



wt^ < 2(4"-^) +«;("-') +4""^)) -4"-^) 



For the remaining inequality, suppose first that wl" = W2 ■ Then 
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If, however, w]^ = w^ , then 

(n) (n) , r„. (n-1) , (n-1) , (n-l)^ (n-l)i (n-1) 

< 2[W2 +'"^3 +W4 ) — Wl — w\ 

, (n) (n) 

< vy -vi ' , 

using the previous inequaUty. This completes the induction step. □ 

The maximum growth rate of the elements at level n of the Apollonian group is also describ- 
able in terms of the joint spectral radius of the generators {5i, 5*2, 5*3, 54} of the Apollonian 
group. 

Definition 7.1. Given a finite set of n x n matrices S = {Mi, . . . , Mg} the joint spectral radius 
cr{T,) is 

(7(1;) :=limsupi max o-(Mj, • • • MjJ^/^ i , 

where a{M) := max{|A| : A eigenvalue of M} is the spectral radius of M. 

The notion of joint spectral radius has appeared in many contexts, including wavelets and 
fractals; see Daubechies and Lagarias ^J for a discussion and references. In general it is hard 
to compute, but here we can obtain an explicit answer. 

Theorem 7.2. The joint spectral radius for the generators S = {Si, S2, S3, S4} of the Apollo- 
nian group is o"(S) = 6^'^ where 



61 = - M + ^/5 + Y 2 + 2\/5 j f« 2.890 . (7.5) 

It is attained by M = S4S3S2S1. 

Proof. Pick a fixed root quadruple with a < 0, say v = (—1, 2, 2, 3), and consider the associated 
packing Vv Lemma 15.51 asserts that 

co||M||f < ||Mv||oo < ci||M||i7, all Me A. (7.6) 

We use the well-known fact that, for any real n x n matrix M, 

a{M) = lim IIM'^II)/'' . (7.7) 

fc— >oo 
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Now ()7.6() gives for any reduced word M = Si^ ■ ■ ■ Si^Si-^ £ A with i^ / ik-i that 

a(M)i/^= lim(||M^v||oo)^ , 
Choosing k = 4n, Theorem 17. II yields 

a{M)^/' < Hm ||r4„^v| 

n— >oo 

Since T^ns = {34,338281)'^^ , this gives 



1 

I 4ns 
00 



1 

I 4ns 



a{M) < hm ||(5453S25i)""v| 
< a{S4S3S2Si)^/^ . 

Choosing M = S4SSS2S1 E A attains equahty (with s = 4), which determines the joint 
spectral radius. A computation reveals that the characteristic polynomial of M = S4S3S2S1 
is X^ - 2X^ - 2^2 - 2X + 1 = which factors as 

{X^ + (-1 + V5)X + 1){X^ - (1 + V5)X + 1) = . 

Its spectral radius is given by (|7.5j) . □ 

The minimal growth rate of any reduced word of length 2n is attained by the word W2n = 
(S^S^)'^. If V = (a, b, c, d) is a root quadruple with a < 0, then the supremum norm 

1 1 W2nv| loo = n{n +l){a + b)-nc+{n- l)d (7.8) 

grows quadratically with n. We omit the easy proof of this fact. 

To conclude this section, we consider the "average value" of ||M^v||oo over all reduced words 
W^ in ^ of length n, which we define to be the median of this distribution. (The elements of the 
distribution are exponentially large, so the median is a more appropriate quantity to consider 
than the mean value.) Let T„ denote the median. We expect that its growth rate should be 
related to the Hausdorff dimension a of the limit set of the Apollonian packing. The results 
of §5 lead to the heuristic that one should expect 

lim llogr„, = i^ . (7.9) 

n^oo n a 

We leave the proof (or disproof) of this as an open problem. 
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8. Open questions 

There remain many open questions concerning integral Apollonian circle packings. We list a 
few of these here. 

(1) In any primitive integral Apollonian packing "P, just four distinct residue classes modulo 

12 can occur as curvature values in V. For example, for V = (0,0,1,1), these values are 
0,1,4,9 (mod 12) while for V = (-1,2,2,3), they are 2,3,6,11 (mod 12). As we noted in 
Section 6, it seems likely that in the packing (—1, 2, 2, 3), all sufficiently large integers congruent 
to 2, 3, 6 and 11 ( mod 12) actually do occur. However, in the packing (0, 0, 1, 1), instead of the 
8 residue classes 0,1,4,9,12,13,16,22 (mod 24) which we might expect to occur, the classes 

13 and 22 are completely missing. Again, computation suggests that only finitely many values 
in the other 6 are missing in (0, 0, 1, 1). Is it true that in any integral Apollonian packing, the 
only congruence restrictions on the curvature values are for the modulus 24? However, in no 
case can we even show that the set of values which do occur has positive upper density. 

(2) Is there a direct way for determining the root quadruple to which a given Descartes 
quadruple belongs? The only way we currently know uses the reduction algorithm described 
in Section 3. 

(3) Concerning root quadruples, what are the asymptotics of the total number of root 
quadruples having Euclidean height below T? 

(4) We have not proved any reasonable lower bound on the number of integers below T 
that occur as curvatures in a fixed integral Apollonian packing. For how large a (3 can one 
prove asymptotically that at least T"~^°^^' integers occur in every such packing? 

(5) All of the preceding questions can also be raised for integral Apollonian packing of 
spheres in 3 dimensions, as discussed in J19j . For example, what are the modular restrictions 
(if any) for the Descartes quintuples (o, 6, c, d, e) occurring in the packing with root quintuple 
(0,0,1,1,1)? 

(6) There exist strongly integral Apollonian packings, in which the circles all have integer 
curvatures and also the curvature x centers of the circles are Gaussian integers, where the circle 
centers are coordinatized as complex numbers. (See J^.) Most questions investigated for 
integral packings can also be asked for strongly integral packings. If we write {x,xX) for a 
circle with curvature x and (complex) center X, then the pairs (x, xX) must also satisfy various 
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modular constraints. For example, modulo 12, the standard integral packing (i.e., starting with 
the circles (—1, 0), (2, 1), (2, —1) has just 20 types of circles, namely, 

{x,xX) = (2,1), (2, 3), (2, 5), (2, 7), (2, 9), (2, 11) 
(3,2i),(3,4i),(3,8i),(3,10i) 
(6, 3 + 4i), (6, 3 + 8i), (6, 9 + 4i), (6, 9 + 8i) 
(11,0), (11,4), (11,8), (ll,6i), (11,4 + 6i), (11,8 + 6i) 

and there are just 120 different four-circle configurations. What are the asymptotics of these 
types and configurations? What is the characterization of the integral (complex) vectors 
{x, xX) that can appear in a given packing? 

We hope to return to some of these issues in a future paper. 
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